An Efficient High-Order Time Integration Method for Spectral-Element Discontinuous Galerkin Simulations in Electromagnetics

نویسندگان

  • Misun Min
  • Paul F. Fischer
چکیده

We investigate efficient algorithms and a practical implementation of an explicittype high-order timestepping method based on Krylov subspace approximations, for possible application to large-scale engineering problems in electromagnetics. We consider a semi-discrete form of the Maxwell’s equations resulting from a high-order spectral-element discontinuous Galerkin discretization in space whose solution can be expressed analytically by a large matrix exponential of dimension κ × κ . We project the matrix exponential into a small Krylov subspace by the Arnoldi process based on the modified Gram–Schmidt algorithm and perform a matrix exponential operation with a much smaller matrix of dimension m ×m (m κ). For computing the matrix exponential, we obtain eigenvalues of the m ×m matrix using available library packages and compute an ordinary exponential function for the eigenvalues. The scheme involves mainly matrix-vector multiplications, and its convergence rate is generally O( tm−1) in time so that it allows taking a larger timestep size as m increases. We demonstrate CPU time reduction compared with results from the five-stage fourth-order Runge–Kutta method for a certain accuracy. We also demonstrate error behaviors for long-time simulations. Case studies are also presented, showing loss of orthogonality that can be recovered by adding a low-cost reorthogonalization technique.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Discontinuous Spectral Element Methods for Time- and Space-Fractional Advection Equations

We develop spectral element methods for a timeand space-fractional advection equation of the form 0Dt u(x, t) + θ 0Dxu(x, t) = f(x, t), of order τ ∈ (0, 1], ν ∈ (0, 1), subject to Dirichlet initial/boundary conditions. We present two spectrally accurate and efficient methods for global discretization of both temporal and spatial terms, instead of employing traditional low-order time-integration...

متن کامل

A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics

In this paper, we discuss the formulation, stability and validation of a high-order non-dissipative discontinuous Galerkin (DG) method for solving Maxwell’s equations on non-conforming simplex meshes. The proposed method combines a centered approximation for the numerical fluxes at inter element boundaries, with either a second-order or a fourth-order leap-frog time integration scheme. Moreover...

متن کامل

GPU Accelerated Discontinuous Galerkin Time Domain Algorithm for Electromagnetic Problems of Electrically Large Objects

In this paper, an efficient time domain simulation algorithm is proposed to analyze the electromagnetic scattering and radiation problems. The algorithm is based on discontinuous Galerkin time domain (DGTD) method and parallelization acceleration technique using the graphics processing units (GPU), which offers the capability for accelerating the computational electromagnetics analyses. The bot...

متن کامل

A Hybridized Crouziex-Raviart Nonconforming Finite Element and Discontinuous Galerkin Method for a Two-Phase Flow in the Porous Media

In this study, we present a numerical solution for the two-phase incompressible flow in the porous media under isothermal condition using a hybrid of the linear lower-order nonconforming finite element and the interior penalty discontinuous Galerkin (DG) method. This hybridization is developed for the first time in the two-phase modeling and considered as the main novelty of this research.The p...

متن کامل

High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere

High-order triangle-based discontinuous Galerkin (DG) methods for hyperbolic equations on a rotating sphere are presented. The DG method can be characterized as the fusion of finite elements with finite volumes. This DG formulation uses high-order Lagrange polynomials on the triangle using nodal sets up to 15th order. The finite element-type area integrals are evaluated using order 2N Gauss cub...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • J. Sci. Comput.

دوره 57  شماره 

صفحات  -

تاریخ انتشار 2013